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摘 要 使 用 多 元 回归 法 进行 调节 效应 分 析 在 社 科 领域 已 常 有 应 用 。 简 述 了 目前 多 元 回归 法 
的 调节 效应 分 析 存 在 的 不 足 , 包括 人 为 变换 检验 模型 、 自 变量 和 调节 变量 区 分 不 足 、 误 差 方 
差 齐 性 的 假设 难以 满足 、 调 节 效 应 量 指标 AR 没有 直接 测量 调节 变量 对 自 变量 与 因 变 量 关系 
的 调节 程度 。 比 较 好 的 方法 是 用 两 水 平 回归 模型 进行 调节 效应 分 析 并 使 用 相应 的 效应 量 指 
标 。 在 介绍 新 方法 和 新 效应 量 后 ， 总 结 出 一 套 调节 效应 的 分 析 流 程 ， 通 过 一 个 例子 来 演示 如 
何 用 Mplus 软 件 进行 两 水 平 回 归 模 型 的 调节 效应 及 其 效应 量 分 析 。 最 后 讨论 了 两 水 平 回归 模 
型 的 调节 效应 分 析 的 发 展 , 包括 稳健 的 调节 效应 分 析 、 潜 变量 的 调节 效应 分 析 、 有 调节 的 中 
介 效 应 分 析 和 有 中 介 的 调节 效应 分 析 等 。 
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1 引言 

调节 变量 (moderator) 在 心理 、 教 育 、 社 会 和 管理 等 研究 中 扮演 着 重要 的 角色 。 如 果 
自 变量 X 与 因 变 量 了 的 关系 受到 第 三 个 变量 Z 的 影响 ， 则 变量 Z 就 是 调节 变量 。 此 时 ，X 
与 了 的 关系 是 调节 变量 Z 的 函数 , 三 者 的 关系 可 以 用 图 1(a) 的 调节 效应 模型 来 表示 。 目前 调 
节 效 应 分 析 常 用 多 元 回归 法 ( 温 忠 己 等 , 2005; 方 杰 等 , 2015)， 但 多 元 回归 法 的 调节 效应 分 
析 及 其 效应 量 指标 AR 仍 存在 不 足 之 处 。Yuan 等 (2014) 提出 两 水 平 回 归 模 型 (two-level 
regression model) 进行 调节 效应 分 析 并 引入 新 的 调节 效应 量 指 标 。 本 文 将 深入 讨论 两 水 平 
归 模 型 的 调节 效应 分 析 及 其 效应 量 。 在 简介 多 元 回归 法 的 调节 效应 分 析 及 其 效应 量 指标 和 
R 后 ， 讨 论 了 多 元 回归 法 和 效应 量 AR? 的 不 足 ， 接 着 详 述 了 如 何 利用 两 水 平 回归 模型 进行 
调节 效应 分 析 并 引入 了 与 此 模型 相应 的 新 效应 量 指标 ， 然 后 总 结 出 一 套 调节 效应 的 分 析 流 
程 。 随后 用 示例 演示 如 何 用 Mplus 软件 进行 两 水 平 回归 模型 的 调节 效应 及 其 效应 量 分 析 , 最 
后 讨论 了 两 水 平 回归 模型 的 调节 效应 分 析 的 发 展 。 
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人) 路 径 图 


图 1 调节 效应 模型 图 (改编 自 方 杰 等 , 2015) 


2 多 元 回归 法 的 调节 效应 分 析 及 其 传统 效应 量 
常见 调节 效应 可 以 通过 下 面 回归 方程 进行 分 析 ( 图 1(5) 是 相应 的 路 径 图 )， 
Y =p, +pX+p,Z+p,XZ +e (1) 
ha B, ARSE, BN“ Z=0 时 自 变 量 X 对 因 变 量 了 的 效应 ，p, 为 当 0 时 调节 变量 Z 对 因 
CO 变量 了 的 效应 ， 假 设 误差 项 e。 服 从 正 态 分 布 。 如 果 系 数 A 的 估计 值 经 统计 检验 显著 ， 则 表 


明 调节 变量 Z 对 自 变 量 X 与 因 变量 了 的 关系 具有 调节 作用 。 调 节 效 应 量 用 AR =R? -R fiih, 


= 


3 方程 C) 去掉 乘积 项 XZ 后 的 测定 系 


N 其 中 Ro 表示 回归 方程 (1) 的 测定 系数 ，Rz 表示 下 


N He GEBZE, 2005, 2020; 方 杰 等 , 2015). 

目前 ， 多 元 回归 法 被 广泛 应 用 于 调节 效应 分 析 。 但 该 方法 仍 存在 诸多 不 足 〈Yuan et al., 
2014)。 第 一 、 多 元 回归 法 将 图 lg) 的 调节 效应 (moderation) 模型 人 为 的 变换 成 1(5) 的 交互 
~ 效应 Gnteraction) 模型 ， 并 通过 最 小 二 乘法 来 估计 此 模型 。 这 种 做 法 其 实 是 检验 Z 与 X 的 交 
= 互 效应 是 否 显著 ， 并 没有 直接 对 图 1(O) 的 调节 效应 模型 进行 检验 。 第 二 、 自 变量 和 调节 变量 
r 的 作用 在 方程 〈1) 中 难以 区 分 。 使 用 多 元 回归 法 ， 调 节 作用 具有 对 称 性 。 也 就 是 说 ， 如 果 
X 对 Y 的 效应 受到 Z 的 调节 《此 时 X 是 自 变量 )， 那 么 Z 对 Y 的 效应 同样 也 受到 X 的 调节 《此 时 
是 调节 变量 )， 且 两 种 调节 效应 的 分 析 结 果 相 同 (Hayes & Montoya, 2017)。 第 三 、 误 差 方差 
齐 性 的 假设 难以 满足 。 多 元 回归 法 假设 误差 方差 齐 性 〈 即 当 自 变量 X 取 不 同 的 值 w 时 ， 误 差 
方差 相等 ) ， 回 归 系 数 显 著 性 的 检验 也 基于 这 一 假设 。 然 而 ， 方 程 (1〉 很 难 满足 误差 方差 
齐 性 的 前 提 假 设 。 因 为 当 自 变量 X 取 茶 值 xz 时 ， 方 程 A) 的 误差 方差 为 


Var(e| X = x,) =Var(Y) -(B, + Bx,)’Var(Z) » 


只 有 当 B，= 0 时 , 误差 的 方差 才 相等 Sx AHR) HALF, 2021) 。Aguinis 等 (1999 ) 
发 现 ， 大约 50% 的 文章 在 使 用 多 元 回归 分 析 时 都 明显 不 满足 误差 方差 齐 性 的 假设 ,此 时 会 导 
MAM p, 的 参数 估计 的 标准 误 有 偏差 ， 从 而 导致 调节 效应 分 析 的 检验 力 Cpower) 下 降 、 效 
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应 量 估计 也 出 现 偏差 〈 刘 红云 等 , 2021) 。 第 四 、 调 节 效应 量 指标 AR 存在 不 足 。AR 并 没 
有 直接 测量 调节 变量 Z 对 X~7 关 系 的 调节 程度 ， 人 入 R? 反 映 的 是 乘积 项 XZ 对 因 变 量 Y 的 方差 的 
额外 解释 力 (Liu & Yuan, in press) 。 温 忠 麟 和 叶 宝 娟 (2014) 认为 人 R? 变 化 要 超过 2% (其 


233%) 才 有 实质 性 意义 , 但 人 RR 的 数值 往往 很 小 (Aguinis et al., 2005; Liu & Yuan, in press) 。 


3 两 水 平 回 归 模 型 的 调节 效应 分 析 及 其 效应 量 


3.1 两 水 平 回归 模型 的 调节 效应 分 析 
两 水 平 回归 模型 的 调节 效应 分 析 ( 图 2) 可 表示 为 (Yuan et al., 2014): 


水 平 1: ¥,=£,+8,X,+e, (2) 
IKF 2: Bo = Yo + YoLi + Ewo (3a) 
pb = Vo t MZ; + en (3b) 


其 中 6 和 忆 是 随机 系数 ， 即 系数 大 小 随 个 体 REMEM, yor Yas Mw 和 yy 为 固定 系 
数 ， 即 系数 大 小 不 随 个 体 i 的 变化 而 变化 。e; Eo e 都 服从 均值 为 0 的 正 态 分 布 。 具 体 


地 ，e ~ N(0,a2)， 加 vol © 7 水 平 1 的 误差 e 和 水 平 2 误差 (eu 和 es ) 相 互 独 


il Go oO; 


AE, 但 eu 和 si 有 可 能 相关 (误差 的 协 方差 Cov(eu,s,) 用 co 表示 )， 自 变量 X 和 调节 变量 Z, 有 


可 能 相关 (图 2 中 用 虚线 表示 )。 将 方程 (3a) 和 (35) 带 入 方程 (2)， 得 : 


Y, =o + Yoli + VX, + HL Xt+En Xt+é0+e, (4) 
如 果 乘 积 项 ZX 的 系数 yi| BEANO, KIZ 对 X, 和 因 变 量 了 关系 的 调节 效应 显著 。 需 要 


说 明 的 是 ， 由 于 方程 (2)-(4) 都 只 有 下 标 i 因此 两 水 平 回 归 模 型 实际 上 是 对 单 水 平 数据 (个 体 
水 平 ) 进 行 调节 分 析 。 也 正 因为 实际 数据 不 是 能 套 的 两 水 平 数据 ， 所 以 两 水 平 回 归 模 型 分 析 


无 法 从 es 中 区 分 出 e ， 也 就 无 法 区 分 cy 和 ac: ， 因 此 令 co =o05 +0. 。 


Level 1 


图 2 两 水 平 回归 模型 的 调节 效应 图 
(改编 自 Yang & Yuan, 2016) 
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两 水 平 回归 模型 用 两 水 平 的 结构 对 单 水 平 数据 进行 调节 效应 分 析 存 在 明显 的 优势 。 第 
一 、 两 水 平 回归 模型 (方程 (2)-(4)， 图 2) 中， 调节 变量 Z, 直接 解释 了 自 变量 X, 和 因 变 量 了 
KHAR H Bo M ) ， 解 决 了 多 元 线性 回归 的 第 一 个 不 足 。 第 二 、 两 水 平 回归 模 型 能 清楚 
地 区 分 自 变 量 X, 和 调节 变量 Z; 的 作用 ， 自 变量 和 在 层 1 方程 中 ， 作 用 是 解释 因 变 量 Y, 的 方 
差 ; 调节 变量 Z, 却 只 能 在 层 2 方程 中 , 作用 是 解释 层 1 截 距 po MRE p DE. HEE X, 
和 调节 变量 Z 是 不 可 互 换 的 ， 它 们 各 自 承担 着 不 同 的 作用 。 解 决 了 多 元 线性 回归 的 第 二 个 
不 足 。 第 三 、 两 水 平 回 归 模 型 不 需要 误差 方差 齐 性 的 假设 , 可 以 在 误差 方差 是 异 方差 的 情况 
下 进行 调节 效应 分 析 。 因 为 ， 两 水 平 回归 模型 (方程 (4)) 的 误差 方差 


Var(e + E£ +€,X,) =X o? +2X,0,, +0, TO =X o +2X,0,,+6,, 


会 随 自 变 量 X 的 变化 而 变化 ， 即 两 水 平 回归 模 型 的 误差 方差 本 就 是 异 方差 的 (Yuan et al., 
2014). Yuan 等 (2014) 的 模拟 研究 也 发 现 ， 当 误差 方差 齐 性 的 假设 不 成 立 ( 即 误差 方差 是 腊 方 
差 ) 时 , 相 比 多 元 回归 法 , 两 水 平 回 归 模 型 的 调节 效应 分 析 可 以 提供 更 准确 的 参数 估计 。Yuan 
等 (2014) 的 模拟 研究 还 发 现 ， 当 误差 方差 齐 性 的 假设 成 立时 ， 两 水 平 回 归 模 型 与 多 元 回归 法 
的 参数 估计 基本 相同 。 这 是 因为 当 s,=0 时， 基于 多 元 回归 法 的 调节 效应 模型 (方程 (1)) 可 看 


WE 


成 是 两 水 平 回归 模型 (方程 (4)) 的 特例 (Liu & Yuan, in press). 
3.2 调节 效应 的 效应 量 
3.2.1 效应 量 指标 : RS 


Yuan “ (2014) 提出 ， 两 水 平 回 归 模 型 的 调节 效应 的 效应 量 (effect size) 指 标 为 


= Var (Yo +7,,Z;) ax f Var(Z,) 


R? ~ ~ => 
Var(B.,) ¥,Var(Z,) +67 


Bil 


tH Var(Z,) 为 调节 变量 Z, PATTIE, Py MST Mo? MHI. R 的 值 即 为 所 的 


Pian 
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方差 中 可 被 调节 变量 2; 解释 的 部 分 ， 而 p, 是 因 变 量 Y, 对 自 变量 X, WIR HVAR 代 


表 了 调节 变量 Z, 对 XX, 与 7 关系 的 调节 程度 。 因 此 ， 用 R3, 作为 调节 效应 的 效应 量 指标 更 符 


合 图 1(O) 中 调节 效应 的 定义 。 大 多 数 情况 下 ，R2 会 比 多 元 回归 法 得 到 的 调节 效应 的 效应 量 


ARK, AKR 表示 的 是 B, 的 方差 被 解释 了 多 少 ， 而 人 R? 表示 的 是 因 变 量 Y 的 方差 被 角 


NI 


RES BZD, ,的 方差 一 般 都 小 于 因 变 量 了 的 方差 (Liu & Yuan, in press) 。 此 外 ， R3, 的 取 
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值 范围 是 [0, 1]。 具 体 地 ，R2 = 1 表示 p 的 变化 完全 是 由 2 的 变化 引起 的 ， 即 完全 调节 ; 0 
< R32 <1, 表示 ,的 变化 部 分 是 由 Z, 的 变化 引起 的 , 即 部 分 调节 ，R3 =0 表示 变量 Z 没 


有 任何 调节 作用 (Yang & Yuan, 2016) 。 综 上 可 知 ， R3, 作为 调节 效应 量 的 指标 相 比 信 R? 


更 为 合适 ， 解 决 了 多 元 回归 法 的 第 四 个 不 足 。 


3.2.2 效应 量 指标 AR MI VR? 系列 


Liu 和 Yuan (in press) 指出 R3, 仍然 存在 不 足 。 R3, 表示 pb 的 方差 中 可 被 调节 变量 2Z, 解 


释 的 部 分 ， 而 不 是 了 的 方差 被 自 变 量 X 解 释 的 部 分 (Var(B,X,) ) 中 可 被 调节 变量 Z, 解 释 的 
部 分 。 为 此 ， 他 们 提出 了 几 种 新 的 调节 效应 量 指标 。 


第 一 种 是 
R 


其 中 Ri 表示 方程 Y=yoo+ yorZ;+ob 的 测定 系数 ， 即 Y 的 方差 能 被 Z; 解 释 的 比例 。R? 表示 


方程 六 =yoo+yorZi+VXioX;+EnXi+o06 的 测定 系数 ， 即 的 方差 能 被 Xx; 和 ZZ; 共同 解释 的 比 


例 。R 表示 方程 (4) 的 测定 系数 ， 即 7 的 方差 能 被 xX;,、Z; 和 ZX; 共 同 解释 的 比例 。AR。 


表示 2Z;X; 独 特 解释 的 7 的 方差 (排除 X; 和 ZZ 的 影响 ) 在 X; 和 2Z,X; 共 同 解释 的 7 的 方差 ( 排 


EZ, 的 影响 ) 中 所 占 的 比例 。AR, 的 取 值 范 围 是 [0, 1]. 


第 三 种 是 所 谓 的 VR? 效 应 量 指标 ， 包 括 
_ Var(y,,Z,X,) _ yrVar(Z,X,) 
Var(BX,) Var(y,+7Z +e)X] > 
yrVar(Z,X,) 
Vary + 7112,)X,] 
y,,Var(Z,X,) 
Varl(y,,Z, +€,)X,]_ 


VR? 


VR = 


和 VR? = 


这 三 个 指标 的 分 子 都 一 样 ,表示 2Z;X; 解 释 的 7 的 方差 ,但 分 母 不 同 。 VR2 和 VR3 的 分 母 是 VR? 


的 分 母 的 一 部 分 , 所 以 VR2 和 VR3 的 值 都 大 于 VR? 。VRI 表示 ZX, 解释 的 了 的 方差 在 与 X; 有 


关 的 了 的 方差 中 所 占 的 比例 。VR2 表示 ZX RERI Y, 的 方差 在 X, 能 解释 的 区 的 方差 中 所 占 


的 比例 。VR3 表示 ZX; 解释 的 7 的 方差 在 Z,X, 解 释 的 了 的 方差 入; 不 能 解释 的 7 的 方差 之 


和 中 所 占 的 比例 ，VR3 类 似 于 方差 分 析 的 效应 量 w; o VR 的 取 值 范围 是 [0, 1], 但 VR? AVR; 
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有 可 能 大 于 1, 此 时 建议 更 换 其 他 效应 量 指 标 或 者 将 自 变量 X 和 调节 变量 Z 都 中 心 化 后 再 进 
行 调节 效应 检验 和 效应 量 计算 。 
3.2.3 效应 量 的 讨论 

效应 量 应 具有 与 测量 单位 无 关 、 不 受 样本 容量 影响 和 


| 


单调 性 的 性 质 ,多 元 回归 法 的 调节 


指标 AR 就 同时 具有 上 述 三 个 性 质 ( 温 忠 鹿 等, 2016) 。AR2 , 具有 与 测量 单位 无 关 和 


sek 
长 
地 


不 受 样本 容量 影响 的 性 质 ， 但 AR;, FE VAS) BEER? -Ri 为 条 件 的 系数 xy 的 单调 函数 (Liu & 


Yuan, in press) 。 当 自 变 量 关 和 调节 变量 Z 都 中 心 化 后 ，VR? VRI 和 VR3 才 具 有 与 测量 单 


位 无 关 和 不 受 样本 容量 影响 的 性 质 ， 但 Liu 和 Yuan 并 没有 讨论 VR? 、VR2 和 VR3 是 否 具有 


单调 性 。 因 为 保持 a 和 4。 不 变 ， axJ(axx+b” 是 x 的 单调 函数 ， PRUAVR? 、VR2 和 VR2 都 有 


单调 性 。Liu 和 Yuan 指出 ， 在 三 个 VR? 效 应 量 指标 中 ，VR? 最 符合 调节 效应 的 定义 。 因 此 ， 


我 们 建议 先 将 自 变量 和 调节 变量 Z 都 中 心 化 ， 再 计算 VR? 效应 量 ， 在 计算 VR? 效应 量 时 ， 


建议 先 计算 VR? ， 并 报告 分 母 的 方差 Var(B,X,) 大 小 ， 便 于 理解 VR? 的 实际 意义 。 然 后 再 根 


据 研 究 者 的 具体 需要 ， 报 告 VR2 MVR o 
4 分 析 流 程 
面 对 一 个 调节 效应 分 析 任 务 , 研究 者 应 当 如 何 进 行 呢 ? 根 据 前 面 的 讨论 , 我 们 总 结 出 一 
套 调 节 效 应 的 分 析 流 程 ( 见 图 3): 
1. 利用 两 水 平 回归 模型 进行 调节 效应 分 析 。 


2. Hre, MATE (3b)) 的 方差 ?是否 显著 。 如 果 是 ， 则 误差 方差 一 定 是 异 方差 的 ， 


报告 两 水 平 回 归 模 型 的 调节 效应 分 析 结 果 和 调节 效应 量 ， 如 果 否 ， 则 表示 不 能 推翻 o? =0 的 


原 假设 ， 还 无 法 得 出 误差 方差 齐 性 ( 即 o?=0) 的 结论 ， 因 此 进入 步骤 3 


3. 利用 多 元 回归 法 进行 调节 效应 分 析 ， 判 断 两 水 平 回归 模型 的 BIC 是 否 小 于 多 元 回归 
法 的 BIC(Bayesian Information Criterion)。 如 果 是 , 则 两 水 平 回 归 模 型 和 数据 拟 合 的 更 好 (Yuan 
et al., 2014)， 报 告 两 水 平 回 归 模 型 的 调节 效应 分 析 结 果 和 调节 效应 量 ， 如 果 否 ， 则 报告 多 元 
回归 法 的 分 析 结 果 和 效应 量 AR“ 。 


两 水 平 回归 模型 的 BIC 


是 否 小 于 多 元 回归 法 ? 


报告 两 水 平 回归 模型 报告 多 元 回归 法 的 分 
的 分 析 结 果 和 效应 量 析 结 果 和 效应 量 AR? 


图 3 调节 分 析 流 程 


5 示例 

接 下 来 用 一 个 实际 例子 演示 如 何 进 行 两 水 平 回归 模型 的 调节 效应 分 析 。 本 例 要 研究 的 是 
冷酷 无 情 特质 ( 思 ) 在 儿童 期 谋 待 (X) 与 青少年 的 网 络 欺负 行为 (7) 之 间 的 调节 作用 , ee 
据 (940 人 ) 均 来 自 方 杰 等 (2020) 的 研究 。 所 有 数据 在 分 析 前 都 已 经 标准 化 。 采 用 Mplus8.2 进 
行 两 水 平 回归 模型 的 调节 效应 分 析 (Mplus 程序 见 附录 )。 需 要 说 明 的 是 ， 第 一 ， 通 过 设置 
每 个 层 2 组 只 有 一 个 层 1 被 试 Mplus 语句 为 CLUSTER = ID), ik Mplus 软件 将 单 水 平 数 
据 当 成 是 两 水 平 数 据 ， 从 而 进行 两 水 平 回归 模型 的 调节 效应 分 析 (Liu et al., 2020)。 第 二 ， 


Mplus 软件 默认 采用 稳健 标准 误 的 极 大 似 然 估计 法 (maximum likelihood estimation with 


= 


robust standard errors ) 进 行 参 数 估 计 ,Liu 等 (2020) 建 议 使 用 贝 叶 斯 法 进行 参数 估计 CMplus8.3 


及 其 以 上 版 本 才能 实现 )。 第 三 ，Mplus 软件 默认 误差 cu 为 0。 因 为 误差 不 相关 (cu =0) 时 ， 


能 得 到 更 准确 的 c2 和 cy 值 (Yuan et al., 2014; Liu et al., 2020, in press; Liu & Yuan, in press). 


两 水 平 回归 模型 的 调节 效应 分 析 结果 见 表 1。 由 于 of 显著 Co? = .29,Z= 3.46)， 因 此 


报告 两 水 平 回 归 模 型 的 调节 分 析 结 果 。 系 数 yi 的 估计 值 显 著 (y， = .11, t=2.89, P= .004), 
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0.085]， 表 明 p, 的 方差 中 有 4% 可 被 调节 变量 ZZ 解释 。 效 应 量 AR2, VR? 、VR2 和 VR2 采用 


R 软件 计算 (Liu & Yuan, in press)。 结 果 表明 ，AR2,= .64， 表 示 ZiX; 独特 解释 的 到 的 方差 


在 X; 和 和 ZX, 共同 解释 的 7 的 方差 中 所 占 的 比例 为 64%; VR? = .043， 表 示 Z,X; 解 释 的 了 的 


方差 在 与 X, 有 关 的 了 的 方差 (.355) 中 所 占 的 比例 为 4.3%; VRI = .229， 表示 ZX 解释 的 也 


的 方差 在 xX; 能 解释 的 7 的 方差 .067) 中 所 占 的 比例 为 22.9%; VR3 = .05， 表 示 Z,X; 解 释 


的 Z 的 方差 在 Z,X; 解 释 的 了 的 方差 和 X; 不 能 解释 的 7 的 方差 之 和 (〈.303) 中 所 占 的 比例 为 
5%, 
表 1 Mplus 软件 的 两 水 平 回归 模型 的 调节 效应 分 析 结 果 


2 2 
Yoo Yoi Yio “1 on oo, R5, 

0 -04 .14 .17 .11 .29 .44 
.04 

SE .03 .03 .03 .04 .08 .08 


Z -1.69 413°" 536°" 289” 3.46” 5.73 


注 :“ 表 示 p<.01; ”表示 p < .001。 


6 两 水 平 回归 模型 的 调节 效应 分 析 的 拓展 
6.1 稳健 的 调节 效应 分 析 

当 数 据 存在 异常 值 或 者 粗 尾 (heavy tails) 时， 两 水 平 回归 模型 的 调节 效应 分 析 结 果 存 
在 偏差 (Yuan et al., 2014; Yang & Yuan, 2016)。 为 了 解决 这 个 问题 ，Yang 和 Yuan (2016) 提 
出 了 两 种 稳健 (robust) 的 两 水 平 回归 模型 的 调节 效应 分 析 方 法 。 一 种 是 基于 t 分 布 权 重 的 
极 大 似 然 估计 法 ， 另 一 种 是 基于 Huber 权重 的 M 估计 法 ， 字 母 “M” 表 示 极 大 似 然 估计 。 
稳健 方法 的 基本 思想 是 ， 为 每 个 数据 分 配 适 当 的 权重 ， 远 离 数据 中 心 的 数据 (如 异常 值 )， 
会 给 予 更 小 的 权重 , 因此 异常 值 对 调节 效应 估计 的 影响 就 更 小 了 ,从 而 得 到 更 准确 的 调节 效 
应 估计 值 。Yang # Yuan (2016) 的 模拟 研究 表明 ， 当 误差 是 非 正 态 (如 粗 尾 ) 时 ， 这 两 种 稳健 
方法 在 调节 效应 分 析 中 的 表现 都 优 于 基于 正 态 分 布 的 极 大 似 然 估计 法 ; 当 误 差 是 正 态 分 布 
时 ， 这 两 种 稳健 方法 在 调节 效应 分 析 中 的 表现 与 极 大 似 然 估计 法 相当 。 

对 示例 进行 稳健 的 两 水 平 调节 效应 分 析 ， 结 果 显 示 ， 两 种 稳健 方法 的 系数 yi 的 估计 值 
都 显著 〈 表 2) ， 表 明 调 节 效 应 都 显著 不 等 于 0， 但 系数 的 估计 值 都 小 于 稳健 标准 误 的 极 大 
似 然 估 计 法 的 分 析 结 果 (yi1=.11)。 


Me 


表 2 稳健 的 两 水 平 回归 模型 的 调节 效应 分 析 结 果 


Huber 权重 1 分 布 权重 
0 > 
6 SE Z 6 SE Z 

Yoo -12 .03 -4.64 -22 .02 -9.81 
Yor 09 .02 3.70 06 .02 2.76 
Yio 19 04 5.53 .17 .03 6.25 
N11 .07 .03 2.72 .05 .02 2.14 
c? .08 .04 1.98 .09 .05 1.86 
al .29 .04 4.25 .17 .03 4.87 
R? 062 025 

Pa i 


6.2 潜 变 量 的 调节 效应 分 析 

显 变量 的 调节 效应 分 析 的 一 个 不 足 是 假设 所 有 变量 的 测量 都 不 存在 测量 误差 , 当 测 量 误 
差 较 大 的 时 候 这 会 造成 调节 效应 的 低估 ， 而 建立 结构 方程 模型 (Structural Equation Model, 
SEM) 进行 调节 效应 分 析 的 最 大 优势 就 在 于 设置 潜 变 量 ， 有 效 控制 测量 误差 ， 能 准确 估计 
调节 效应 值 ， 是 比较 好 的 方法 。 


具体 地 ， 假 设 潜 变量 由 p, 个 指标 ,x%,…, x 组 成 ， 潜 变量 5 由 p, MEM ccs, 


HR WEE ny, 由 p, 个 指标 y,y,,…,y,, 组 成 , Liu 等 (2020) 指出 , 两 水 平 回归 模型 的 SEM 


调节 效应 分 析 就 是 将 方程 (2) 和 (3) 中 的 显 变量 了 XM DIARRE, Eu FUE, » 


同时 将 图 2 中 的 长 方形 框 (表示 显 变 量 ) 换 做 椭圆 形 图 框 (表示 潜 变 量 ) 即 可 。 调 节 效 应 量 指标 


R3, 只 需要 将 方程 (8) 的 Var (Z,) ACA Var (é) 即 可 。 


Liu 等 〈2020) 的 模拟 研究 发 现 ， 当 误差 方差 是 异 方差 时 ， 相 比 单 水 平 的 SEM 调节 效 
应 分 析 ( 包 括 乘积 指标 法 和 潜 调 节 结 构 方 程 法 ) 而 言 ， 两 水 平 回归 模型 的 SEM 调节 效应 分 忆 
提供 了 更 准确 的 调节 效应 估计 、 且 可 靠 区 间 和 覆盖 率 和 第 I 类 错误 率 的 表现 都 更 优 (覆盖 率 更 
接近 95%， 第 I 类 错误 率 更 接近 0.05); 当 误 差 方 差 齐 性 时 ， 两 水 平 回归 模型 的 SEM 调节 效 
应 分 析 和 单 水 平 的 分 析 结果 相当 。Liu 等 提供 了 Mplus 程序 进行 两 水 平 回归 模型 的 SEM. 调 


a 


节 效 应 分 析 。 值 得 注意 的 是 ， 目 前 还 不 清楚 AR。。 和 VR? 系列 指标 是 否 适用 于 两 水 平 回 归 模 


型 的 SEM 调节 效应 分 析 (Liu & Yuan, in press) 。 


6.3 基于 两 水 平 回 归 模型 的 有 调节 的 中 介 效 应 分 析 

当前 ， 调 节 和 中 介 的 整合 已 经 成 为 社 科 领域 发 展 的 趋势 。Liu 等 (in press) 提出 了 基于 
两 水 平 回归 模型 的 有 调节 的 中 介 效 应 的 分 析 方 法 。 例如 ,基于 两 水 平 回归 模型 的 后 半 中 介 路 
径 (MY) 被 调节 的 中 介 效 应 分 析 可 表示 为 如 下 方程 : 


水 平 1: Y,=8,+bM,+cX, +e, (5) 
水 平 2: Bu =VotVoL: + En (6a) 

b, =o t+ V,Z, + Eq (6b) 
水 平 ]: M,= £6, +aX,+e,, (7) 


WR ya 显著 不 为 0( 即 yiia 的 置信 区 间 不 包含 0)， 就 表示 中 介 效 应 被 调节 了 。Liu 等 Gn 
press) 将 被 调节 的 中 介 效 应 的 效应 量 定义 为 中 介 效 应 ab, = alo + yZ, + En) 的 方差 中 , BEB Z, 


解释 的 比例 ， 即 [y a’ Var(Z )]/Varlab,) =[y,a’°Var(Z)I/Ly,a’Var(Z,) + a’Var(e,, | > FII, 读者 很 


容易 得 到 ， 基 于 两 水 平 回归 模型 的 前 半 中 介 路 径 (XM) 被 调节 的 中 介 效 应 分 析 可 表示 为 
如 下 方程 : 


水 平 1: M,= By, + 4,X, + ey; (8) 
水 平 2: Bui = Vw HZ: + En (9a) 

a, =o t+ M,Z, + Eq (9b) 
水 平 1: Y =p, +bM,+c'X, +e, (10) 


如 果 yib 显 善 不 为 0( 即 yb 的 置信 区 间 不 包含 0)， 就 表示 中 介 效 应 被 调节 了 。 被 调节 的 中 
介 效 应 的 效应 量 定义 为 中 介 效 应 ab= (o +y Z +6 ,)b 的 方差 中 ， 能 被 2, 解 释 的 比例 ， 即 


[y2b’Var(Z,)]/Var(a,b) =[y}b’°Var(Z )I/ly7b’Var(Z)+b°Var(e,)] o Liu 等 提供 了 Mplus 程序 进行 


两 水 平 回归 模型 的 有 调节 的 中 介 效 应 及 其 效应 量 分 析 。 

Liu 等 (in press) 的 模拟 研究 发 现 ， 当 误差 方差 是 异 方差 时 ， 相 比 单 水 平 的 有 调节 的 中 介 
模型 而 言 , 两 水 平 回归 模型 的 有 调节 的 中 介 效 应 分 析 能 提供 更 准确 的 条 件 中 介 效 应 估计 、 上 有 
FY Se DX TA) ain SL 类 错误 率 的 表现 都 更 优 ( 履 盖 率 更 接近 95%， 第 I 类 错误 率 更 接近 
0.05); 当 误 差 方 差 齐 性 时 ， 两 水 平 回归 模型 的 有 调节 的 中 介 效 应 分 析 和 单 水 平 的 分 析 结 果 
相当 。 
6.4 基于 两 水 平 回归 模型 的 有 中 介 的 调节 效应 分 析 

除了 基于 两 水 平 回归 模 型 的 有 调节 的 中 介 效 应 分 析 ， 刘 红云 等 (2021) 还 提出 了 基于 两 
水 平 回 归 模 型 的 有 中 介 的 调节 效应 ( 即 调节 变量 Z 对 XY 关系 的 调节 作用 通过 中 介 变 量 MM 
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a 


来 传递 ) 的 分 析 方 法 。 基 于 两 水 平 回归 模型 的 有 中 介 的 调节 效应 分 析 《〈 图 4) 可 表示 为 如 下 


方程 : 
水 平 1: 了 = +cX +e， (11) 
水 平 2: Peet Va + 72 + En (12a) 
C= Ni + IM, +y: + En (12b) 
水 平 2: M, =p, +aZ,+eé,, (13) 
水 平 2 
2 a reece eee eee eee 
> 水 平 1 
一 图 4 两 水 平 回归 模型 的 有 中 介 的 调节 分 析 
NI (改编 自 刘 红 云 等 , 2021) 
与 两 水 平 回归 模型 的 调节 效应 分 析 (方程 (2) 和 (3)) 相 比 ， 两 水 平 回归 模型 的 有 中 介 的 


A 调节 效应 分 析 仅 增加 了 方程 (13) 并 且 在 方程 (12) 中 增加 了 中 介 变 量 M, 而 已 。 将 方程 (12a) 


= Al (12b) 带 入 方程 (11)， 得 : 
€ Y, = Vo yao M; +yoZ: + MX i +y M, X; +722 X; tE X, + Ey; (14) 


其 中 ，s, = s+e,， 因 为 在 单 水 平 数据 中 (方程 44) 中 仅 有 下 标 门 ，sio Me, 无 法 区 分 开 的 。 将 


方程 (13) 带 入 方程 (14)， 整 理 得 
Y, = Vo + Yobu +o + yo)Z; + [po Fd + (Vna ty)2 + (Ve + En)]X, + ie Ey 
其 中 ，yiia+yis 表 示 Z 对 XY 关系 的 调节 作用 ， 其 中 yy, 是 直接 调节 作用 ，yija 是 Z 通 过 


M 对 X— 了 关系 的 调节 效应 ， 即 有 中 介 的 调节 作用 。 如 果 yiia 显著 不 为 0( 即 yiia 的 置信 区 


间 不 包含 0)， 就 表示 Z 对 X-> 了 关系 的 调节 作用 被 M 中 介 了 。 刘 红云 等 (2021) 将 有 中 介 
的 调节 效应 的 效应 量 定义 为 自 变 量 X 对 因 变 量 了 的 效应 c; 的 方差 Var(c) 中 ， 能 被 有 中 介 的 


Val 13 380M ya 解释 的 比例 ， 即 [zazVar(Z )]/Varc) 。 刘 红云 等 提供 了 Mplus 程序 进行 两 水 平 


回归 模型 的 有 中 介 的 调节 效应 及 其 效应 量 分 析 。 
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刘 红 云 等 (2021) 的 模拟 研究 发 现 ， 当 误差 方差 是 异 方差 ( 即 s, 关 0) 时 ， 相 比 单 水 平 的 有 
中 介 的 调节 模型 (去 掉 方程 (12) 的 &, ) 而 言 ， 两 水 平 回归 模型 的 有 中 介 的 调节 效应 分 析 能 提供 
更 准确 的 参数 估计 和 效应 量 估 计 、 第 I 类 错误 率 更 接近 0.05， 并 且 随 着 &, 的 增 大 ， 两 水 平 
回归 模型 的 优势 随 之 增加 ; 当 误 差 方 差 齐 性 ( 即 es, =0) 时 ， 两 水 平 回归 模型 的 有 中 介 的 调节 效 


应 分 析 和 单 水 平 的 分 析 结 果 相 当 。 


7 结语 

针对 多 元 回归 法 的 调节 效应 分 析 及 其 效应 量 AR 的 不 足 (Yuan et al., 2014; Liu & Yuan, 
in press)， 较 好 的 改进 方法 是 使 用 两 水 平 回归 模型 进行 调节 效应 分 析 并 使 用 与 此 模型 相应 的 
新 效应 量 指标 。 我 们 总 结 出 一 个 调节 效应 分 析 流 程 , 并 通过 一 个 实例 演示 了 如 何 进行 两 水 平 
回归 模型 的 调节 效应 分 析 并 计算 其 效应 量 , 随后 还 讨论 了 两 水 平 回 归 模 型 的 调节 效应 分 析 的 
发 展 。 但 本 文 仍然 存在 一 些 不 足 ， 尚 需 进 一 步 深入 讨论 和 拓展 。 例 如， 本 文 仅 涉及 一 个 自 变 
一 个 调节 变量 和 一 个 因 变 量 的 调节 效应 分 析 ， 实际 上 ， 本 文 介绍 的 方法 同样 适用 于 多 个 
自 变量 、 多 个 调节 变量 的 情况 。Yuan 等 人 (2014) 用 两 水 平 回 归 模 型 分 析 了 两 个 自 变 量 、 三 个 
调节 变量 的 调节 效应 。Liu 等 (in press) 用 两 水 平 回归 模型 分 析 了 两 个 调节 变量 分 别 调节 前 、 
后 路 径 的 有 调节 的 中 介 效 应 ,方法 的 进步 给 研究 者 提供 了 一 个 深入 理解 和 应 用 两 水 平 回 归 模 
型 进行 调节 效应 分 析 的 机 会 相信 随 着 研究 的 深入 , 会 不 断 增加 我 们 对 两 水 平 回 归 模 型 的 调 
节 效 应 分 析 问 题 的 理解 。 
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Abstract: In recent years, multiple regression has been widely used in social sciences to analyze 
the moderating effect. However, this practice was found to have at least four weaknesses. First, the 
concept of moderation is artificially treated as interaction. Second, the role of the predictor X is 
confounded with that of the moderator Z. Third, the assumption of homoscedasticity of error 
variances across different values of X and Z is often violated by data in social and behavioral 
sciences. Fourth, AR? do not directly measure the effect of moderation as conceptually defined. 
Compared to multiple regression, the two-level regression model has many advantages in the 
analysis of moderating effect. After introducing the moderation analysis method based the 
two-level regression model and the corresponding effect size, we propose a procedure to analyze 
the moderating effect based on the two-level regression model. We exemplify how to conduct the 
proposed procedure by using Mplus. Directions for future study on two-level regression models 
are discussed at the end of the paper, including robust method of two-Level regression models; 
two-level moderated latent variable model, two-level moderated mediation model, and two-level 
mediated moderation model. 


Key words: moderating effect, two-level regression model, multiple regression, effect size 


附录 : 两 水 平 回 归 模 型 的 调节 效应 分 析 的 Mplus 程 序 

DATA: FILE = 2.txt; 

VARIABLE: NAMES ARE ID x zy; 1!2.fxt 的 变量 名 
USEVARIABLES = x z y xz; 
CLUSTER =ID; DÆ 
WITHIN = x zy xz; ! 水 平 1 的 变量 
DEFINE: xz=x*z; ! 产 生 调节 项 xz 

ANALYSIS: TYPE IS TWOLEVEL RANDOM; 

ALGORITHM = INTEGRATION; 


MODEL: 

JoWITHIN% 
b| Y ON X; ! 设 置 随机 斜率 
YONZ lgam01 值 

XZ(r11);”!gam11 值 

X with Z; 
X with Y; 

%BETWEEN% 
[b]; ! gam10 值 
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b(sig2_eil); ! sigmal 的 平方 
MODEL CONSTRAINT: NEW(var_z R2); 
var_z=1; !Z 由 于 标准 化 ， 方 差 为 1 
R2=(r11*r11*var_z)/((r11*r11*var_z)+sig2_eil1); ! 求 效应 量 RS 
OUTPUT: CINTERVAL; 
TE: 本 程序 默认 使 用 MLR 估 计 ， 如 果 使 用 贝 叶 斯 估计 ， 则 需要 将 ANALYSIS 命 令 改 为 : 
ESTIMATOR IS BAYES; 
POINT =MEAN;  ! 后 验 均 值 作为 参数 的 点 估计 值 


